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Abstract 

The KARMEN experiment uses the reaction ^^C ( fe , e~ ) ^^Ng.s. to measure 
the energy distribution of Ve emitted in muon decay at rest — > e+ + z^e + '^/x- 
The z^e analog ujl of the famous Michel parameter p has been derived 
from a maximum-likelihood analysis of events near the kinematic end point, 
Enr^.^ = 52.8 MeV. The resuh, ujl= (2.7t|^ ± 3.1) x lO^^, is in good agree- 
ment with the standard model prediction lol = 0. We deduce a 90% 
confidence upper limit of ujl < 0.113, which corresponds to a limit of 
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I 9rl + ^S'rl I— 0-'''^ interference term between scalar and tensor 

coupling constants. 

Experimental results from nuclear j3 decay and muon decay form the basis of the V-A 
hypothesis, which is an essential feature of the standard model (SM) of electroweak interac- 
tions. The rate of muon decay, the purely leptonic process ji^ e'^ + Ue + t^^, has been used 
to determine the universal Fermi coupling constant Gp- Precise measurements of the shape 
of the e+ energy spectrum, the decay asymmetry between the /i"*" spin and momentum, 
and the polarisation vector of the have led to bounds on the scalar, vector and tensor 
coupling constants, which form the Lorentz structure of the charged weak interaction. These 
results combined with the inverse process + e" underpin the SM assumption 

of lepton number conservation, the V-A interaction and universality [jl| . All experiments up 
to now support the V-A structure of the weak interaction; however, substantial non- (V-A) 
components are not ruled out. 

Complementary to these experiments, which are all based on observation of the charged 
leptons only, the Karlsruhe Rutherford Medium Energy Neutrino experiment (KARMEN) 
determines the energy spectrum of the v^. emitted in the decay yU^ ^ -|- z^e + of un- 
polarized muons to draw conclusions on the Lorentz structure. In the well-known case of 
e"*" spectroscopy, it is the Michel parameter p which governs the shape of the e+ energy 
spectrum. In an analogous way, the shape of the Vf. energy spectrum is determined by the 
parameter tu^, which also depends on vector, scalar, and tensor components of the weak 
interaction, but in a different combination. In the SM all non- (V-A) components vanish, 
and is predicted to be 0. Thus an upper hmit on cul derived from the analysis of the 
energy spectrum provides new limits on nonstandard couplings. 

All features of muon decay are most generally described by a local, derivative- free, lepton- 
number-conserving, four-lepton point interaction with the matrix element given by 
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M=4^Gp 9l{ee\T'^\il^e)n){il^,)m\TM. (l) 

e,iJ, — R,L 

The index 7 labels the type of interaction F (S = 4-scalar, V = 4-vector, T = 4-tensor) and 
the indices e and /i indicate the chirality (L = left-, R = right-handed) of electron and muon 
spinors, respectively. In this representation the chirality of the neutrino n or m is fixed to 
be equal to that of the associated charged lepton for the V interaction, but opposite for the 
S and T interactions. As Gp sets the absolute strength of the interaction, the ten coupling 
constants g]^ are dimensionless complex quantities normalized by 

3i^Lr+3i^7Lr+ E i\\9^\' + \9r/) = I (2) 

€,ij,=R,L 

with g]ij^ = gj^i^ = 0. In the SM, muon decay is a pure V interaction mediated between 
left-handed particles, so all coupling constants vanish except g^i = 1- Although this rep- 
resentation is elegant from the theoretical point of view, the individual coupling constants 
cannot be determined directly by experiment. However, the measurable parameters (p, t], 
ul, etc.) are expressable as positive semidefinite bilinear combinations of g^^ from which 
upper or lower limits for the coupling constants can be derived. 

The possibility of measuring ul with the KARMEN experiment was first pointed out by 
Fetscher [|]. More recently Greub et al. [@ have calculated the spectrum of left-handed 
including radiative corrections and effects of finite lepton masses. Taking significant terms 
only, the spectrum dNi/dx can be described by 

^ = {G'o(x) + G,{x) + u,G,{x)] (3) 

where is the muon mass, x = 2Ej,/'m^ is the reduced neutrino energy, and denotes 
the probability of emission of a left-handed Ug. The function Go{x) describes the pure V-A 
interaction, Gi{x) takes into account radiative corrections, and ujlG2{x) includes the effect 
of scalar and tensor components according to 

, ^3 \9L\' + M9U' + \9'rl + ^91l\' ... 

' \9k? + \9L? + M9L? + M9lR? + l2\gl,r ^ ^ 



The calculated Ue energy spectra for different values of ujl are shown in Fig. ||(a). Mo- 
mentum conservation in the decay fixes the emission direction of v^. near the kinematic end 
point to be opposite to that of the positron and the z/^. Together with angular momentum 
conservation this implies suppression of emission of left-handed Ve in the case of vector cou- 
pling, while all other couplings enhance the decay rate at the end point. The total decay 
rate, and therefore the integral neutrino flux, is unchanged by nonstandard interactions. 

The KARMEN experiment uses the pulsed spallation neutron facility ISIS at the Ruther- 
ford Appleton Laboratory to investigate neutrinos from decay. The 800 MeV proton beam 
from ISIS is stopped in a Ta-D20 target producing neutrons and pions. All charged pions 
are stopped inside the target within 10~^° s, the vr" being absorbed by the heavy target 
material while the vr"^/ /i+ decay chain 7r+ — * yU+ + z/^, /i+ e+ + z/g + produces an intense 
burst of z/^, z/g, and z/^, emitted isotropically with equal intensity. Since both vr'^ and /i"*" 
decay at rest, the energy spectra of the neutrinos are well defined. The 7r+ decay produces 
monoenergetic with Ey^ = 29.8 MeV; the u^. and z/^ from the fi~^ decay have continuous 
energy distributions up to -Emax = 52.8 MeV. The time structure of ISIS — two 100 ns wide 
proton bunches 324 ns apart and recurring at 50 Hz — determines the production time of the 
different u flavors: the short vr^ lifetime {t^^ = 26 ns) leads to two z/^ pulses within the first 
500 ns after beam-on target. These pulses are well separated in time from the production of 
z/g and z/^, which follow with the much longer lifetime of the fi^ {'^n^ ~ This leads 

to a suppression factor of about 10^ for cosmic-ray background. 

The neutrinos are detected in a segmented 56 ton liquid scintillation calorimeter consist- 
ing of 512 optical modules, each with a length of 3.53 m and a cross section of 18 x 18 cm^ [Q. 
The detector is an almost completely (96%) active calorimeter optimized for the measure- 



ment of electrons around 30 MeV and achieves resolutions of a{E)/E = 11.5%/y -E(MeV) 
for energy, and cr{X) ^ 7 cm for position measurement. A 7000 ton shielding steel block- 
house together with two layers of active veto counters suppresses beam-correlated spallation 
neutrons and cosmic-ray muons. 

The signature that unambiguously identifies a z/g is a delayed coincidence consisting of 
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an electron from the charged current reaction ^^C {ue ,e~ ) ^^Ng.s. in the time window of 
z/g production followed by a positron from the subsequent /3 decay of ^^Ng.s. (r = 15.9 ms) 
at the same location in the detector. Each event fully contained within the central detector 
with time 0.6-9.6 fis after beam-on target and energy 10-36 MeV is identified as electron, 
provided it is followed by a positron event within 0.5-36 ms with energy 3.5-16.5 MeV. 
We demand the sequence to be detected in the same or adjacent module within a distance 
AX < 35 cm along the module axis. Cuts used to reduce cosmic background are the same as 
used in previous data evaluations p . In data accumulated between June 1990 and December 
1995 — corresponding to 9122 C of protons or 2.51 x 10^^ fi^ decays in the ISIS target — 
we find 513 e~ /e^ sequences. Subtracting 13.3 ± 0.8 background events and taking into 
account an overall detection efficiency e = 32.8%, the fiux-averaged cross section is 

(cr)^^p = (9.4 ± 0.4(stat) ± 0.8(sys)) x 10"^^ cml (5) 

This is in good agreement with different theoretical calculations of (cr)^j^ in the range of 
(9. 1-9.4) X 10-^2 i^g]^ 

As the recoil energy of the ^^N nucleus is negligible, the energy is determined from 
the measurement of the electron energy E^. via the kinematic relation E^ = E^ + Q, where 
Q = 17.3 MeV is the Q value of the detection reaction. The energy dependence of the cross 
section is dominated by the phase-space factor [E^ — QY. Therefore, a low rate of additional 
z/g at the kinematic end point -Emax = 52.8 MeV due to nonstandard couplings is translated 
to the observation of a significantly higher rate of electrons and thus to a distortion of the 
visible energy spectrum of Fig. |l](b). 

The KARMEN calorimeter allows a precise measurement of the energy E^ [see Fig. ||(c)]. 
The energy spectrum of v^. from /i"^ decay is then determined in two steps. First, we apply the 
procedure of regularized unfolding described by Blobel to derive the true electron energy. 
This method takes into account the detector response and minimizes inherent instabilities 
(oscillating solutions) by demanding a priori a certain degree of smoothness of the true 
electron distribution depending on statistical accuracy. The Ve energy distribution is then 
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calculated from the number of primary electrons, within a given interval Ai? from the 
unfolding procedure, divided by the corresponding mean cross section. This yields a Vf. 
energy spectrum with seven data points as shown in Fig. ^ and compared with the V-A 
expectation. This represents the first measurement of the neutrino energy spectrum from 
muon decay in addition to the well-known e"*" spectrum. 

Because of the strong energy dependence of the detection cross section the most detailed 
information on ujj^ and Q\ is obtained from the experimental electron spectrum of Fig. |l](c). 
The analysis is done by two independent methods: (1) the investigation of the measured 
decay rate on the basis of the flux-averaged cross section, and (2) the analysis of the spectral 
shape with a maximum likelihood (ML) method. 

As can be seen from Fig. [l|(b), cul > would result in additional ^^C ( z/g , e~ ) ^^Ng.s. 
events; on the other hand, Q\ < 1 would reduce the number of events [see Eq. (H)]. In 
order to find allowed regions in the Q^-ul parameter space, we compared measured and 
expected flux averaged cross sections. As theoretical cross section (o")^^ with a realistic 
estimate of the systematic error we use (o")^^ = (9.2 ± 0.5) x 10"^^ cm^. The experimental 
cross section is taken from Eq. (|^) with statistical and systematic error added quadratically. 
The probability distribution of the ratio 

(cr) 9 4 ± 9 

RiQl, ^l) = ^ = QUI + S.uj,) = ^--^ (6) 
Wth 9.2 ±0.5 

incorporates a flux decrease by right handed through Q1 as well as an increase by nonzero 
ul values; S is the ratio of additional events in case of ujl = 1 relative to the expectation 
for ul = 0. We have sampled the probability density function of the ratio R from Gaussian 
distributions of (o")p^p and {cr)^_^ for 3 different energy ranges: (a) the range 10-36 MeV with 
the highest statistical accuracy, but only moderate sensitivity S = 0.81, (b) the range 28- 
36 MeV, where with S = 3.48 we are very sensitive to ul, and (c) the range 10-22.5 MeV, 
where the expected event number is almost independent of ul [S = 0.002, see Fig. |l|(b)]. 
From range (c) we deduce a lower limit > 0.796. The shaded parameter space shown in 
Fig. ^ combines regions excluded at 90% confidence level of all 3 energy ranges. From inverse 
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muon decay experiments it is known that Q1 > 0.92 [^,0. Including this information in 
our analysis of range (b) restricts the allowed area and sets a 90% confidence upper limit 
ul < 0.12. 

In the second method we determine a;^ by analyzing the shape of the visible electron 
spectrum independent of Q^. In order to increase the energy resolution and to reduce the 
background level we applied more stringent cuts on the electron position along the module 
axis |X| < 150 cm and on the electron time 0.6-7.2 /is. These cuts reduce the background 
to only 6.0 events in a sample of 441 events, thus nearly doubling the signal-to-background 
ratio. 

The theoretical Ug energy spectrum of Eq. |^ was converted into a visible electron spectrum 
using the energy-dependent ct{E^) taken from folded with the detector response by a MC 
calculation. The ML procedure was carried out on an event-by-event basis for several fit 
intervals all of which gave results compatible with ujl = within a la-error. The net result 
is 

= (2.7l|^(stat) ± 3.1(sys)) x 10"^. (7) 

Including the systematic error (energy shift of 0.25 MeV or 0.7% scaling error) we find, with 
the most conservative Bayesian approach, a 90% confidence upper limit ul < 0.113. This 
excludes the region above the horizontal line in Fig. |^. Combining Eq. (H) and Eq. (Q) the 
following relation between the shape parameter ul and nonstandard couphngs is |ll0|,p!2| 
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\9rl + ^91l\<^Y^l. (8) 

The limit on ul thus results in an upper limit of | gf^^ + 2 g]^i^ \ < 0.78 for the interference 
term of scalar and tensor amplitudes. 

In conclusion, the KARMEN experiment finds no evidence for nonstandard coupling 
constants in /x^ decay at rest, either by a determination of the absolute Ug flux or by analysis 
of the spectral shape. Our analysis excludes most of the Q^-wl parameter space and yields 
for the first time an upper limit on the neutrino Michel parameter ui. 



During 1996 the experiment was upgraded by an additional active veto counter in order 



to increase the sensitivity of the search for neutrino oscillations in the channel u^^Ue |]T3 
Since 1997 KARMEN has been taking data again. Up to the end of 1999 we expect about 
400 further charged current events, which will reduce the statistical error by about a factor 
of 1.4. Considering also a reduction of the systematic error, this may result in a limit 
competitive with the present best limit | gf^^j^ + 2 g]^j^ \ < 0.45 deduced from measurements 
of the positron polarization ||T0[| . 
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FIG. 1. Influence of different values of ujl = 0.0, 0.15, 0.3 on (a) the Ue energy spec- 
trum in yU"*" decay and on (b) the visible electron energy spectrum measured with the reaction 
""^■^C ( i^e , e~ ) ""^^Ng.s. (c) Experimental electron energy distribution together with MC expectation 
(solid line) and the subtracted background (shaded). 
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FIG. 2. Energy spectrum of v^. from decay determined by an unfolding method compared 
with the standard model expectation (solid line). 
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FIG. 3. The Q^-ul parameter space: The shaded regions are excluded at 90% confidence from 
the different analyses of the absolute flux in several energy ranges. The horizontal line is the result 
of the spectral shape analysis ujl < 0.113 at 90% confidence. The vertical line is the current best 
limit Ql > 0.92. 
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